Imaging nano-scale electronic inhomogeneity in lightly doped Mott insulator 

Ca2_a;Na^Cu02Cl2 
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The spatial variation of electronic states was imaged in the lightly doped Mott insulator 
Ca2-a:NaxCu02Cl2 using scanning tunneling microscopy / spectroscopy (STM/STS). We observed 
nano-scale domains with a high local density of states within an insulating background. The ob- 
served domains have a characteristic length scale of 2 nm ('-^4-5a, a : lattice constant) with preferred 
orientations along the tetragonal [100] direction. We argue that such spatially inhomogeneous elec- 
tronic states are inherent to slightly doped Mott insulators and play an important role for the 
insulator to metal transition. 

PACS numbers: 74.72.-h, 68.37.Ef, 74.50.+r 



Holes and electrons doped into antiferromagnetic Mott 
insulators can bring about many exotic phenomena such 
as high-temperature superconductivity (HTS) in copper 
oxides and colossal magnetoresistance (CMR) in man- 
ganese oxides [l|. When Mott insulators are doped 
lightly, a conflict arises between the itinerancy of the 
doped charge carriers and the energy gain of the magnetic 
background, which can cause inherently spatially inho- 
mogeneous electronic states, such as stripes 12*, |3, 0, 0, • 
It has been suggested that suclian_ inhomogeneity is a key 
ingredient of HTS and CMR 



Indeed, spatially inhomogeneous superconductivity 
has been observed in Bi2Sr2CaCu20y by scanning tun- 
neling microscopy / spectroscopy (STM/STS) |llialil|. 
It is not yet clear however whether or not the observed 
inhomogeneity is electronic in origin or related to the lo- 
cal off-stoichiometry. This is, in our opinion, because 
the materials used in these STM studies were heavily 
doped superconductors, whereas purely the electronic 
phase separation is anticipated to occur near the doping- 
induced insulator to metal transition (IMT) found in the 
lightly doped region. It is therefore highly desirable to 
extend such real-space imaging into very lightly doped 
region. Bi2Sr2CaCu20y, however, is not chemically sta- 
ble when lightly doped and is therefore an unsuitable 
choice of material to extend STM/STS studies towards 
the IMT. 

An alternative candidate with which to explore the 
real space physics of lightly doped Mott insulators is 
the cupric oxychloride Ca2-xNaajCu02Cl2. The par- 
ent compound Ca2Cu02Cl2 is an antiferromagnetic Mott 
insulator. By substituting Na for Ca, the system is 
doped with holes and eventually develops into a super- 
conductor at X ~ 0.08 jl^i- Moreover, single crystals 
of Ca2-a;Naa;Cu02Cl2 are very easy to cleave, making 
them ideal systems for STM/STS studies. As shown 
in Fig. QJ), Ca2-a^Naa^Cu02Cl2 has a crystal structure 



closely related to those of La2-a;Sra;Cu04 but maintains 
undistorted tetragonal Cu02 layers down to low tem- 
peratures. Such a simple structure with small num- 
ber of constitution elements minimize the possibility of 
the local off-stoichiometry other than doping. In this 
Letter, we report atomic resolution STM/STS results 
on Ca2-a?Naa;Cu02Cl2 with various hole concentrations 
across a doping-induced IMT. A variation of electronic 
states, in the form of nano-scale domains with spatially 
preferred orientations, was clearly resolved. We suggest 
that the observed electronic inhomogeneity is inherent to 
slightly doped Mott insulators. 

Sizeable single crystals of Ca2-a;Naa;Cu02Cl2 were 
grown using a novel high-pressure flux technique over a 
wide concentration range across the IMT boundary [I3 • 
We have carried out STM/STS measurements on single 
crystals from x = 0.06 (insulator) to x = 0.10 (heav- 
ily underdoped superconductor with Tc = 13 K) under 
ultrahigh vacuum (UHV) conditions of ~ 10~^ Pa. The 
samples were cleaved along (001) planes in UHV at about 
80 K and immediately transferred to the STM head and 
cooled to 4.7-7.0 K, where all reported data were taken. 
The exposed surface layer is believed to be the highly 
ionic Ca-Cl layer, since bonding to the adjacent Ca-Cl 
layer is weakest at this interface. 

A typical constant-current STM image of x = 0.08 
crystal, which is in the critical vicinity of the IMT bound- 
ary, is shown in Fig. ^. Regular square atomic lattices 
are clearly observed, and the lattice spacing is estimated 
to be 0.39 nm, in good agreement with the in-plane lat- 
tice constant of the tetragonal unit cell 13]. Superim- 
posed on the regular atomic lattices, we observe nano- 
scale corrugations, which can be seen as the bright (high) 
and dark (low) regions in Fig. ^. The corrugations are 
not spatially regular but rather show a preferred orienta- 
tion with respect to a specific crystallographic axis. As 
shown in Fig. ^, the domains in many cases are aligned 
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FIG. 1: A 15 nm X 15 nm STM images of (001) cleav- 
age planes of Cai.92Nao.o8Cu02Cl2 single crystals taken at a 
sample bias (Vs) of 200 mV, a tunneling current (It) of 10 pA 
and at 7.0 K. (b) Crystal structure of Ca2-ccNaxCu02Cl2. (c) 
Auto-correlation patterns of (a). 



along [100] or [010], and form river-like pathways. Some 
images appear more patchy, but again the patches are 
aligned along [100] or [010]. Auto-correlation analyses 
shown in Fig. ^ clearly and quantitatively demonstrate 
an enhanced spatial correlation along the [100] and [010] 
directions. Around the origin, the correlation shows a 
rapid decay on the length scale of ~2 nm, which corre- 
sponds to the width of the rivers or the size of the patches. 
A weak tail then follows, representing the river-like spa- 
tial correlation. 

We ascribe the observed corrugations to a spatial vari- 
ation of electronic states, rather than a topographic ef- 
fect. It is emphasized here that the height change, de- 
tected by constant-current imaging, represents not only 
the 2;-axis coordinates of the surface atoms but also the 
spatial variations of the local density of states (LDOS). 
If there exists a region with a lower LDOS than its sur- 
roundings, the STM tip will come close to the surface 
in order to keep the tunneling current constant. During 
our STM measurements, we found that the magnitude 
of the corrugations strongly depends on the bias volt- 
age. This implies that the observed corrugations can- 
not be understood simply in terms of a topographic ef- 
fect, but rather by a spatial variation of the energy de- 
pendence of the LDOS. Such spatial variation in LDOS 
was indeed confirmed by the STS measurements shown 
in Fig. I2t- By comparing Fig. ^ and Fig. [^b, it is 
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FIG. 2: (a) dl/dV spectra measured along the line shown 
in (b). Each spectrum is shifted vertically dependent on the 
position. For clarity, the spectra taken at apparently lower 
(higher) positions are colored blue (red), (b) STM image 
with the trajectory along which the spectra in (a) were mea- 
sured. The colored- line code corresponds to the positions of 
the spectra individually-colored shown in (a). 



clear that there is a strong correlation between the shape 
of the differential conductance {dl/dV) spectra, which 
is proportional to the LDOS, and the apparent height. 
This demonstrates that the electronic states of under- 
doped Ca2-a:Naa^Cu02Cl2 are spatially inhomogeneous 
and suggests the formation of distinct domains in the 
local electronic structure. It may be interesting to infer 
that the length scale of the domains observed here is com- 
parable to the size of superconducting islands observed 
in heavily doped Bi2Sr2CaCu20^ !9!,lialil|- 

The dl/dV spectra in Fig. EK do not represent the 
magnitude of the LDOS, due to the change in the tip- 
sample distance in constant current imaging. However, 
if the tunneling barrier is uniform, we can easily esti- 
mate a correction associated with the change in the tip- 
sample distance z, since dl/dV = C ex.p{—2hzz)LDOS^ 
where C is a constant and hz is related to a tunneling 
barrier \\J^. The tunneling barrier height was measured 
along the line shown in Fig. Eb and was found to be 
uniformly 3.3 ± 0.2 eV. We thus performed a numerical 
correction of the tip-sample distance and estimated the 
relative magnitude of the LDOS at each position. Fig. 
indicates the comparison of the relative LDOS averaged 
over the bright and dark regions in the inset of Fig. 01 
It is evident from this plot that LDOS is substantially 
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FIG. 3: The averaged LDOS from the high (red, sohd hne) 
and the low (blue, broken line) domains are shown for a 
X = 0.08 sample. The LDOS was estimated from the dl/dV 
spectrum by calibrating the tip-sample distance. The verti- 
cal ticks located it 10 meV indicate the knees in the LDOS 
spectra. The inset shows the topographic image taken at 
Vs — 200 mV and It ^ 100 pA with markers where the spectra 
were taken (red cross=bright / blue diamond=dark). Colors 
of markers correspond to the colors of spectra in the main 
panel. 



suppressed in the dark regions up to a few hundred meV. 
Therefore, the dark region is substantially less metallic 
(more insulating) than the bright region. We ascribe 
the V-shaped behavior of LDOS to the large pseudo- 
gap (~0. 1-0.2 eV) universally observed in underdoped 
cuprates by angle-resolved photoemission spectroscopy 
(ARPES) near (7r,0). Since this large pseudogap be- 
comes larger upon approaching the undoped Mott insula- 
tor 15], the gap enhancement in the dark regions implies 
the suppression of hole density. Besides the large pseu- 
dogap, a knee-like structure is found around ±10 meV at 
the bottom of the V-shape. We suspect that this might 
originate from superconductivity. 

As the Ca and CI atoms are strongly ionic, their elec- 
tronic states are energetically far away from the Fermi 
level. We therefore believe that the majority of the tun- 
neling current observed originates from the uppermost 
Cu02 plane just below the exposed Ca-Cl layer, and that 
the electronic inhomogeneity takes place within the Cu02 
plane. It might be argued that the observed phenomena 
are surface rather than bulk effects. However, this system 
is extremely two-dimensional in nature and the coupling 
between the adjacent Cu02 planes is negligibly small. It 
is therefore likely that the electronic states of the topmost 



Cu02 layer, sandwiched by ionic Ca-Cl layers, resemble 
those in the bulk. Recent transport measurements in- 
dicate that the anisotropy between the in-plane and the 
out-of-plane resistivities reaches 10^ even in the super- 
conducting X = 0.10 composition [l^. This is one of the 
largest anisotropic resistivity ratios ever observed in the 
cuprate superconductors. 

We consider two scenarios for the origin of observed 
electronic inhomogeneity within the Cu02 plane. One 
is a poorly screened impurity potential. The dopant 
Na+, substituted for Ca^+, should modify the electro- 
static potential of the neighboring Cu02 planes, which 
can give rise to inhomogeneous distribution of carriers. 
In ordinary metals, however, the electrostatic potential 
is screened by itinerant carriers over a subatomic spac- 
ing p^ and can not simply lead to a nanoscale inho- 
mogeneity. This situation may be altered significantly 
in lightly doped cuprates. A close proximity to the 
Mott insulator can bring about a nonlinear screening 
effect \\A [l^l or a reduction of the itinerant coherent 
states l2i|. As a result, effective charge screening length 
can be anomalously long and a nano-scale electronic in- 
homogeneity may emerge. 

Another possible scenario for the observed inhomo- 
geneity may be the electronic phase separation. Elec- 
tronic phase separation arising from competing ordered 
states is at the heart of the physics of CMR phenomena 
in the manganites [^ |^ , where nano-scale phase separa- 
tion between a ferromagnetic metal and a charge ordered 
insulator has been observed by STM/STS HI HI H^. 
Another strong case for electronic phase separation may 
be made for the lightly doped cuprates, and this has been 
discussed to be intimately connected to the mechanism 
of high temperature superconductivity 6] . If this is the 
case, naively, two distinct phases should be observed and 
the ratio of the two should change upon doping. 

To examine this point, we performed STM/STS on 
the samples with different dopings and found system- 
atic evolution of both STM images and dl/dV spectra 
as shown in Fig. ^ It is clear as shown in Fig. ^-c that 
the nano-scale mixture of high and low LDOS regions 
is observed regardless of the ground state, insulator or 
superconductor. The high LDOS (bright) area increases 
with increasing density of holes, and with further doping, 
it is likely that the electronic inhomogeneity will eventu- 
ally fade out. The system appears to change from an 
insulator to a metal once the metallic domains become 
connected over the entire surface. Although these appear 
to be consistent with the phase separation, the evolution 
of the electronic states within the two regions can not 
be understood in its naivest form. As shown in Fig. ^ 
and Fig. ^, while the shape of dl/dV spectra in the 
dark regions are nearly identical, those in the bright re- 
gions show a substantial change as a function of doping, 
which is not consistent with "naive" phase separation 
picture. Even if a kind of electronic phase separation is 
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FIG. 4: Doping dependence of STM images and dl /dV spec- 
tra of Ca2-xNaxCu02Cl2. The topographic images (a), (b), 
and (c) were from x = 0.06 (insulator), x = 0.08 (right at 
IMT), and x = 0.10 (superconductor) samples, respectively. 

(a) was taken at V; = 300 mV, h ^ 100 pA, and T = 5.0 K. 

(b) and (c) were taken at V; = 200 mV, h ^ 100 pA, and T 
= 4.7 K. The images from x — 0.08 (b) was taken on different 
surface from that of Fig. |3 demonstrating the reproducibility 
of main features, (d) and (e) are averaged differential con- 
ductance spectra collected from the bright and dark regions, 
respectively. The spectra are shifted vertically for clarity and 
the horizontal bars indicate zero. The vertical tick marks in- 
dicate the energies of the shoulder structures in the spectra. 



responsible for the observed phenomena, additional in- 
gredient, for example, poorly-screened impurity poten- 
tial, should be invoked. In any event, we believe that the 
observed nano-scale electronic inhomogeneity represents 
the physics of the lightly doped Mott insulator. 

In conclusion, STM/STS measurements on under- 
doped Ca2-a:Naa^Cu02Cl2 rcvcalcd pronounced elec- 
tronic inhomogeneity near the IMT. The observed fea- 
tures of inhomogeneity, i.e. a long characteristic length 
scale of ~2 nm and the evolution of the LDOS as a 
function of doping, can not be ascribed to the conven- 
tional charge screening nor to the naive electronic phase 
separation. We infer that competing orders inherent to 



doped Mott insulators plays an important role for the 
observed electronic inhomogeneity. Regardless the ori- 
gin, real space physics of lightly doped Mott insulator 
is indispensable in understanding the insulator to metal 
transition in a real system. 

We thank N. E. Hussey, N. Shannon, and K. M. Shen 
for critical reading of manuscript. This work was partly 
supported by a Grant-in- Aid for Scientific Research from 
the MEXT, Japan. 

note added in proof: Spatial modulations of tunnel- 
ing conductance with checkerboard-like patterns have re- 
cently observed in the same compound. This is to be 
published elsewhere. 
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